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ABSTRACT 

Using the APP coronagraph of VLT/NACO we searched for planetary mass companions around 
HD115892 and HD172555 in the thermal infrared at 4 jum. Both objects harbor unusually luminous 
debris disks for their age and it has been suggested that small dust grains were produced recently in 
transient events (e.g., a collision) in these systems. Such a collision of planetesimals or protoplanets 
could have been dynamically triggered by yet unseen companions. We did not detect any companions 
in our images but derived the following detection limits: For both objects we would have detected 
companions with apparent magnitudes between ~13. 2-14.1 mag at angular separations between 0.4- 
1.0" at the 5-cr level. For HD115892 we were sensitive to companions with 12.1 mag even at 0.3". 
Using theoretical models these magnitudes are converted into mass limits. For HD115892 we would 
have detected objects with 10-15 Mj up at angular separations between 0.4-1.0" (7-18 AU). At 0.3" 
(—5.5 AU) the detection limit was >25 Mj up . For HD172555 we reached detection limits between 
2-3 M Jup at separations between 0.5-1.0" (15-29 AU). At 0.4" (-11 AU) the detection limit was 
>4 Mj up . Despite the non-detections our data demonstrate the unprecedented contrast performance 
of NACO/APP in the thermal infrared at very small inner working angles and we show that our 
observations are mostly background limited at separations >0.5". 

Subject headings: stars: formation — planets and satellites: formation — planets and satellites: 
detection — stars: individual (HD172555, HD115892) — planetary systems 



1. INTRODUCTION 

While most dedicated surveys to directly image ex- 
trasolar planets around nearby stars yielded null results 
(e.g., Chauvin et al. 2010; Heinze et al. 2010; Lafreniere 
et al. 2007a; Kasper et al. 2007), some remarkable ex- 
ceptions were discovered in the last years: The HR8799 
planetary system (Marois et al. 2008, 2010), Fomalhaut 
b (Kalas et al. 2008), (3 Pictoris b (Lagrange et al. 
2009a; Lagrange et al. 2010), and 1RXS J1609-2105 b 
(Lafreniere et al. 2008, 2010). The host stars of the first 
three systems are all A-type stars and they harbor both 
massive planets and debris disks. For Fomalhaut and j3 
Pictoris dynamical interactions between the exoplanets 
and the disks led to observable signatures: an offset be- 
tween the debris disk center and the star (Fomalhaut) or 
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1 Based on observations collected at the European Organisation 
for Astronomical Research in the Southern Hemisphere, Chile, un- 
der program number 060.A-9800(J). 



disk warps (/3 Pictoris). Thus, although there seems to 
be no direct correlation between the existence of both a 
debris disks and (an) exoplanet(s) (e.g., Apai et al. 2008; 
Moro-Martfn et al. 2007), specific properties of the debris 
disk can hint toward the existence of low-mass compan- 
ions. 

Here, we report on the search for low-mass companions 
around HD115892 and HD172555 using direct imaging. 
Both objects are also young, nearby, A-type stars (see, 
Table 1) that are surrounded by debris disks (Oudmai- 
jer et al. 1992; Moor et al. 2006; Su et al. 2006; Morales 
et al. 2009). These objects are, however, particularly in- 
teresting because their disks have a very high fractional 
luminosity and appear too luminous for their age (Moor 
et al. 2006; Wyatt et al. 2007). Since steady state evo- 
lutionary models of debris disks predict a much lower 
dust luminosity it was suggested that either the plan- 
etesimals in these systems have unusual properties or the 
observed, luminous dust could be transient (Wyatt et al. 
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2007). Interestingly, Lisse et al. (2009) found evidence 
for Silica dust and SiO gas in the mid-infrared spectrum 
of HD172555 which could be indicative of a high-velocity 
collision of protoplanets or planetesimals. Such a colli- 
sion could have been triggered via dynamical interactions 
or gravitational stirring by a so far unseen companion, 
so we sought to search for it directly. 

We used the Apodizing Phase Plate (APP) corona- 
graph installed at VLT/NACO (Kenworthy et al. 2010; 
Girard et al. 2010). The APP is designed to work in the 
3-5 /ttm wavelength range where it enhances the contrast 
between ^2-7 X/D on one side of the PSF (Figure 1; sec 
also, Kenworthy et al. 2007; Codona et al. 2006). This 
inner working angle (IWA) corresponds to projected sep- 
arations of ~5-25 AU around our targets, comparable to 
the giant planets' orbits in our own Solar System. The 
debris disks are located within the inner 10 AU around 
each object (Table 1). 

2. OBSERVATIONS AND DATA REDUCTION 

The data were obtained on 2010-04-04 during the com- 
missioning of the APP with the high-resolution AO- 
camera NACO (Lenzen et al. 2003; Rousset et al. 2003) 
mounted on ESO's VLT UT4. Using the same observ- 
ing setup already used to image the exoplanet j3 Pic- 
toris b (Quanz et al. 2010), we chose the L27 camera (~ 
27.15 mas pixcU 1 ) with the visible wavefront sensor. All 
images were taken in the NB4.05 filter (A c = 4.05 ^m, 
AA = 0.02 /im) and in pupil stabilized mode. These are 
the first data sets that combine the APP with Angular 
Differential Imaging (ADI) (Marois et al. 2006) (see be- 
low). We used the "cube mode" readout where all image 
frames, i.e., each single exposure, are saved individually. 
To ensure that no frames were lost we only read out 
a 512x512 pixels sub-array of the detector. The effec- 
tive field-of-view (FoV) using the APP in this sub-array 
mode is restricted to the uppermost 512x90 pixels (i.e., 
roughly 13.9" x 2.4") as the APP introduces a vertical 
shift of the image along the detector's y-axis 2 . For both 
sources, several data cubes were taken each at a slightly 
different dither position following a 3 point dither pat- 
tern along the x-direction of the detector's effective FoV. 
Halfway through the observations the camera was ro- 
tated by 180° so that the high contrast side of the APP 
covered both hemispheres. Due to an error during the 
rotation of the camera, we did not cover the full 360° 
around the targets (see, section 3). In total we obtained 
54 and 36 data cubes for HD115892 in hemisphere 1 and 
2, respectively, and 18 and 24 data cubes for HD1 72555. 
Each cube consists of 200 individual image frames, i.e., 
exposures. Table 2 summarizes the observations and also 
the observing conditions. To enhance the signal-to-noise 
(S/N) of potential companions, we chose to saturate the 
core of the stellar PSFs, but we note that the APP re- 
duces the peak flux in the PSF core by roughly 40% 
(Kenworthy et al. 2010). 

For the photometric calibration we also obtained un- 
saturated images of both targets. We used the same 
observing strategy but decreased the detector integra- 
tion time (DIT) to 0.0558 s for HD115892 and 0.2 s for 
HD172555 and took only 6 data cubes, each consisting 
of 200 exposures, for each calibration data set. In Fig- 

2 See, NACO User Manual v. 88 page 53. 



ure 1 we show the median combined image of 1 cube for 
HD172555. 

The general data reduction approach (bad pixel correc- 
tion, sky subtraction) is described in Quanz et al. (2010). 
This time, however, since we had sufficient field rotation 
(>10°) during our observations we used the LOCI al- 
gorithm (Lafreniere et al. 2007b) to subtract the stellar 
PSF of our images. LOCI creates a reference PSF for 
each image from a linear combination of all other images 
observed at a different parallactic angles. The coefficients 
of the combination are optimized inside different sub- 
sections of the image independently so that the residual 
noise within each subsection is minimized. We refer the 
reader to the original paper for a more detailed descrip- 
tion of LOCI. By scaling and inserting the PSF of the 
unsaturated images as fake planets with known bright- 
ness in the raw frames and retrieving them in the final 
image we did a small parameter study to optimize the 
LOCI parameters for our purposes. The best results in 
terms of planet contrast and S/N 3 were achieved with the 
following LOCI parameters which we used for the final 
analyses: FWHM=4.5 px, ^=0.75, dr=3, N A =2QQ. We 
note that we used LOCI on each individual image frame 
and not on stacked images as the former approach pro- 
vided better detection performances. The data for each 
hemisphere was reduced separately. In a last step, we 
cut vertically through the center of each PSF-subtracted 
frame and saved only the high-contrast, i.e., left-hand, 
side of the PSF. For both hemispheres these images were 
then rotated to the same field orientation and averaged 
to create our final image. No additional filtering was 
applied to our data. 

3. RESULTS & ANALYSIS 

In Figure 2 we show the final PSF-subtracted images 
for HD115892 and HD172555 in the left column. Due 
to the rotation error described above we lack data in a 
wedge in the North-East quadrant of both objects. On 
the opposite side of the central star, however, a wedge 
of the same size was covered during the observations of 
both hemispheres. The right column of Figure 2 shows 
the number of frames that were eventually combined for 
a given position for the final image. Our final images 
probe regions as close as ~0.3" (5.5 AU) and ^0.4" (11 
AU) around HD1 15892 and HD172555, respectively, and 
out to 1" as the maximum distance around each target. 
Since objects tend to slightly drift across the detector in 
NACO's pupil stabilized mode, we could only combine 
the maximum number of frames in the innermost ~0.8" 
(see, right column Figure 2). 

We did not detect any faint companions in our final 
images. However, by inserting and retrieving fake plan- 
ets we can determine the sensitivity of our observations 
and put constraints on the maximum brightness of po- 
tential non-detected companions. Since LOCI can sig- 
nificantly reduce the flux of any detected point source, 
detection limits need to be based on the retrieval of fake 
companions and cannot be derived solely from the noise 
in the final image. A complicating factor in our case 
is the inhomogeneous sensitivity across the final images. 
We decided to put fake planets in regions where only 

3 A description how we define planet contrast and S/N is given 
in section 3. 
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data from one hemisphere is combined and not in the 
overlapping regions. This approach is representative for 
a "typical" APP observing run and the results are rep- 
resentative for typical APP detection limits. We used 
the data of hemisphere 1 for HD115892 and of hemi- 
sphere 2 for HD172555. For HD115892, 10681 and 9615 
frames were combined for separations <0.8" and >0.8", 
respectively. For HD172555 we could combine 4697 and 
3046 frames in these regions. As fake planets we used 
for each target the median-combined PSF of an unsatu- 
rated data set and scaled it to different contrast ratios 
based on the average count rate of the unsaturated im- 
ages and the difference in exposure time between the un- 
saturated and the saturated images. These fake planets 
with known brightness were then inserted in the individ- 
ual sky-subtracted raw frames at different radii taking 
into account the field rotation that occurred between the 
exposures. Finally, we repeated the data reduction de- 
scribed above and determined the S/N of fake planets 
that we recovered in the final image. We did aperture 
photometry on the recovered planets and compared it to 
the standard deviation of background pixels in an annu- 
lus centered on the central star. This annulus had the 
same radial distance as the planet and a width twice as 
wide as the aperture radius. We excluded those regions 
in the annulus where fewer frames were combined than 
at the position of the planet, and we excluded the region 
around the planet itself (i.e., 3 FWHM centered on the 
planet) as LOCI can create artificial 'holes' left and right 
of a detected point source. The S/N of the fake planet 
can then be expressed as 

S/N = F pl /(a ■ (1) 

with F p i being the flux of the planet, a the standard 
deviation of the pixels in the annulus (both measured in 
'count rate') and r ap the aperture radius. We inserted 
fake planets with a contrast between 9 and 11 mag in 
the HD 115892 data and between 8 and 9 mag in the 
HD 172555 data and computed the S/N for two aperture 
sizes (2 and 3 pixels radius). The final 5a contrast limit 
for a given separation was then derived by averaging the 
S/N in both apertures, taking those fake planets where 
the averaged S/N was the lowest but >5 and extrapo- 
lating the contrast of the planet to a value that would 
correspond to a 5a detection. We emphasize that we did 
not apply any sort of filtering or background smoothing 
to our data which makes our final S/N estimates rather 
conservative. Also the optimized extraction of a PSF 
template could lead to the robust detection of fainter 
companions. 

In Figure 3 we show the final 5a detection limits 
for both objects between 0.3-1.0". Overplotted are de- 
tectable mass limits for a given contrast and the age of 
the star (Table 1). These mass limits are derived from the 
DUSTY and COND evolutionary models (Chabrier et al. 
2000; Baraffe et al. 2003). We use the COND models for 
objects with effective temperatures below —1700 K and 
the DUSTY models for hotter objects. For the 350 Myr 
old object HD115892 our data reach a contrast between 
—10.5-11.3 mag at angular separations between 0.4-1.0" 
(7-18 AU). This contrast corresponds to detectable mass 
limit between 10-15 M Jup . At 0.3" (-5.5 AU) the con- 
trast is —9.4 mag and we are still sensitive to objects 



with masses >25 M Jup . For the 12 Myr HD172555 sys- 
tem the contrast is —9.2-9.8 mag at separations between 
0.5-1.0" (15-29 AU) which corresponds to mass limits 
of 2-3 M,j up . At 0.4" (-11 AU) the achieved contrast is 
—8.9 mag and we are still sensitive to objects with >4 
Mjup. Due to the smaller field rotation for this object 
we can not probe IWA <0.3". 

Both our datasets have comparable total integration 
times and factoring in the apparent brightness of the 
stars both curves are comparable in terms of detectable 
brightness for potential companions. In addition, both 
contrast curves are relatively flat for separations >0.5". 
This suggests that the APP achieves close to background 
limited performance for these separations. We computed 
the expected background limit for the HD115892 data set 
based on the sky noise in individual frames far away from 
the star. The dashed line in the left panel in Figure 3 
shows the result and confirms that our data are indeed 
(mostly) limited by the background and not by the con- 
trast. Due to the lack of appropriate dark frames we 
could not repeat this exercise for HD 172555. 

Given the non-homogeneous data coverage in azimuth 
the detection limits derived above vary between differ- 
ent positions around each object. To estimate the global 
detection limits we only consider those regions where we 
have combined at least half as many frames as for the 
analysis. For HD115892 we then have to exclude an az- 
imuthal wedge between — 22°-134° (East of North), and 
for HD172555 it's regions between -17°-107° (see also 
right column Figure 2). In all the other parts of the im- 
ages the detection limits shown in Figure 3 apply with 
a significance of >3.5cr with the lowest significance ap- 
plying only in very small wedges directly adjacent to the 
excluded parts. 

4. DISCUSSION 

Lagrange et al. (2009b) used radial velocity to search 
for planetary mass companions to both of our targets. 
Although they did not find any they could put some 
constraints on the occurrence of massive planets in short 
period orbits. For HD115892 they could exclude objects 
more massive than 1.7, 3.8 and 100.0 Mj up in 3, 10, 
and 100 day orbits, respectively, with >99% confidence. 
These orbital periods correspond to semi-major axes of 
—0.06, 0.13, and 0.62 AU assuming circular orbits. For 
HD 172555 the same confidence level was achieved for ob- 
jects more massive than 11.3 Mj up in a 3 or 10 day orbit 
(i.e., -0.05 AU and -0.11 AU), respectively. For a 100 
day orbit (i.e., at 0.53 AU) the detection limit was 33 
M Jup . In addition, Biller et al. (2007) used NACO in 
spectral differential imaging (SDI) mode to search for 
low-mass companion around HD172555. While our data 
are more sensitive in the innermost 1" (i.e., for separa- 
tions <30 AU), their data cover regions out to 2" and 
they could have detected objects with masses >5 Mj up 
for separation between —30-60 AU. 

In combination with the other studies our data put 
stringent constraints on the existence of giant exoplan- 
ets around HD 115892 and HD 172555 interior and exte- 
rior to the debris disks. Although our data probe re- 
gions very close to the assumed location of the debris 
disks (4-6 AU), our hypothesis that dynamical interac- 
tions between a planet and the debris disk could have 
led to a recent collision of planetesimal lacks direct ob- 
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servational support. We note, however, that planets (or 
planetary systems) with masses below our detection lim- 
its are certainly able to dynamically shape debris disks 
and influence their evolution (see, e.g., Raymond et al. 
2011). 

Our data demonstrate that the APP opens up a new 
parameter space for direct imaging of exoplanets by 
pushing the background limit significantly closer to the 
star. A comparison to surveys carried out in the H band 
shows that Chauvin et al. (2010) and Lafreniere et al. 
(2007a) typically reached a contrast of —10 mag and —9.5 
mag at a separation of 0.5", respectively. For HD115892 
our contrast is >11 mag at the same separation. A sim- 
ilar contrast has been reported by the NICI campaign 
at the Gemini observatory (Chun et al. 2008) also op- 
crating in the H band, but a more direct comparison of 
the contrast performance is limited due to different in- 
tegration times and different target stars with different 
brightnesses. However, since planetary mass objects ap- 
pear red in the infrared, NACO/APP has an advantage 
when it comes to the final detectable mass limits, because 
it works in the L band and not in the H band. 

5. SUMMARY & CONCLUSIONS 

We presented the first observations combining NACO's 
Apodizing Phase Plate coronagraph with Angular Dif- 
ferential Imaging to search for faint companions to the 
young debris disk host stars HD 115892 and HD 172555 in 
the NB4.05 filter. Our conclusions are as follows: 

• We did not detect any point sources but achieved 
the following detection limits: For HD115892 we 
could have detected objects with a contrast of 
—10.5-11.3 mag (corresponding to 10-15 Mj up ) at 
angular separations between 0.4-1.0" (7-18 AU). 
At 0.3" (—5.5 AU) the detection limit was a con- 
trast of -9.4 mag (>25 M Jup ). For HD172555 we 
reached a contrast of —9.2-9.8 mag (2-3 Mj up ) at 
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separations between 0.5-1.0" (15-29 AU). At 0.4" 
(—11 AU) the detection limit was a contrast of 8.9 
mag (—4 Mj up ). These limits are >3.5ct limits. We 
do not have data in an azimuthal wedge between 
-22°-134° (East of North) for HD115892 and be- 
tween ~17°-107° for HD172555. 

• While current/previous high-contrast imaging 
campaigns carried out in the H band are contrast 
limited at small IWA, our data are mostly back- 
ground limited for separations >0.5" in the thermal 
infrared. 

• Taking advantage of the red H-L color of planetary 
mass objects, NACO/APP is capable of detecting 
cooler planets (i.e., lower mass or older planets) 
compared to observations in the near-infrared for a 
given contrast. 

NACO/APP is currently a superior combination to 
search for planets at unprecedented small IWA in par- 
ticular around bright targets. And even when the next 
generation high-contrast imaging instruments such as 
SPHERE and GPI come online, with its unique L-band 
capabilities NACO/APP can help to characterize at least 
a certain subset of the exoplanets these instrument will 
find in the near-infrared. 
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Fig. 1. — The PSF of HD172555 observed with NACO/APP. The image shows the median combination of one cube of unsaturated 
exposures. In the left-hand side the APP effectively suppresses the diffraction rings between ~2-7 A/D. 
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Fig. 2. — Final PSF subtracted images of HD115892 and HD172555 (left column) and corresponding images showing the number of 
frames that were combined for a given position (right column). All images have a linear scale. The pixel units of the PSF subtracted 
images is count rate and the stretch ranges from -5cr to +5<r where a denotes the standard deviation of the counts in the background. We 
computed a in those regions in the image where more than the mean number of frames where combined (i.e., 7710 frames for HD115892 
and 3494 frames for HD172555). Due to an error during the rotation of the camera we lack data for a wedge in the North-East quadrant 
around both targets (North is always up and East to the left). 
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Fig. 3. — 5a detection limits (squares) given as magnitude contrast in the NB4.05 filter as a function of radial separation from the host 
star for HD115892 (left panel) and HD172555 (right panel). The lower x-axes show the radial separation in arcsec while the upper ones 
depict the projected separation in AU. Due to the larger field rotation in the dataset for HDlf5892 we can probe inner working angles 
as small as 0.3". Overplotted are the expected contrast for planets with different masses (dotted lines). For HDf 15892 we also plot the 
measured background limit for our observations (dashed line). See text for more details. 



TABLE 1 

Basic properties of target stars based on Su 
et al. (2006) and wyatt et al. (2007) and 

REFERENCES therein. 



Parameter 



HD115892 



HD172555 



RA (J2000) 
DEC (J2000) 
Apparent mag 
Spectral type 
Mass 
Age 

Distance 

Debris disk radius 



in L 



13' l 20 m 35 s .82 

-36°42'44".26 

2.68 mag a 

A2V 

2.5 M 

350 Myr 

18 pc 

6 AU 



18' l 45 m 26 s .90 
-64°52'16".53 
4.28 mag b 
A5IV/V C 
2.0 Mq 
12 Myr 
29 pc 
4 AU 



a Morcl & Magncnat (1978) 

b Based on K-band magnitude from Cutri ct al. 
(2003) and a K-L color of 0.02 mag for an A5V star (see: 
http: / / www.jach.hawaii.edu /UKIRT / astronomy /utils / temp.html) . 

c Gray et al. (2006) found a spectral type of A7V. 



TABLE 2 

Summary of deep imaging observations in pupil tracking mode. 



Parameter 


HD115892 


HD115892 


HD172555 


HD172555 




1st hemisphere 


2nd hemisphere 


1st hemisphere 


2nd hemisphere 


UT start 


02h:48m: 19.28s 


05h:03m:52.30s 


06h:44m:42.80s 


08h:26m:32.36s 


UT end 


04h:43m:13.37s 


06h:16m:41.86s 


08h:05m:19.54s 


10h:llm:23.80s 


NDIT x DLT a 


200 X 0.5 s 


200 X 0.5 s 


200 X 1.2 s 


200 X 1.2 s 


NINT b 


54 


36 


18 


24 


Parallactic angle start 


-78.99° 


-10.04° 


-71.76° 


-43.78° 


Parallactic angle end 


-28.52° 


53.63° 


-48.88° 


-7.81° 


Airmass 


1.19. . . 1.03 


1.02. . . 1.05 


1.71. . . 1.46 


1.42. ..1.31 


Typical DIMM seeing [A=500 nm] 


0.6". . .0.8" 


0.5". . .0.7" 


0.5"... 0.6" 


0.5". . .0.9" 


(-EC)mcan / (-EC*) m i n / {-EC) max C 


46.2 / 8.5 / 58.6 % 


49.0 / 14.6 / 64.1 % 


36.97 / 23.3 / 48.8 % 


47.1 / 25.2 / 58.9 % 


(ro)mcan / (T()) min / (r ( )) max d 


9.0 / 5.2 / 13.1 ms 


9.3 / 5.6 / 13.0 ms 
78.95°... 141. 73° 


4.1 / 2.9 / 5.3 ms 


6.0 / 3.7 / 8.5 ms 
45.21°... 79. 76° 


PA c 

r "camera 


-170.00°... -121.00° 


-162.77°. . .-141.06° 



a NDIT = Number of detector integration times (i.e., number of individual frames); DIT = Detector integration time (i.e., single 
frame exposure time). 
b NINT = Number of data cubes. 

c Average, minimum and maximum value of the coherent energy of the PSF in data cubes. Calculated by the Real Time 
Computer of the AO system. 

d Average, minimum and maximum value of the coherence time of the atmosphere in data cube. Calculated by the Real Time 
Computer of the AO system. 
°Position angle of camera adaptor at the beginning of exposure. 



